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Efficient feedback controllers for continuous-time quantum error correction 
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We present an efficient approach to continuous-time quantum error correction that extends the 
low- dimensional quantum filtering methodology developed by van Handel and Mabuchi [quant- 
ph/05 11221 (2005)] to include error recovery operations in the form of real-time quantum feedback. 
We expect this paradigm to be useful for systems in which error recovery operations cannot be ap- 
plied instantaneously. While we could not find an exact low-dimensional filter that combined both 
continuous syndrome measurement and a feedback Hamiltonian appropriate for error recovery, we 
developed an approximate reduced-dimensional model to do so. Simulations of the five-qubit code 
subjected to the symmetric depolarizing channel suggests that error correction based on our ap- 
proximate filter performs essentially identically to correction based on an exact quantum dynamical 
model. 
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I. INTRODUCTION 

Quantum error correction is inherently a feedback pro- 
cess where the error syndrome of encoded qubits is mea- 
sured and used to apply conditional recovery operations 
[1]. Most formulations of quantum error correction treat 
this feedback process as a sequence of discrete steps. Syn- 
drome measurements and recovery operations are per- 
formed periodically, separated by a time-interval chosen 
small enough to avoid excessive accumulation of errors 
but still comparable to the time required to implement 
quantum logic gates [1, 2]. There is, however, mounting 
evidence from the field of real-time quantum feedback 
control [3, 4, 5, 6] that continuous observation processes 
offer new, sometimes technologically advantageous, op- 
portunities for quantum information processing. 

Toward this end, Ahn, Doherty and Landahl (ADL) [7] 
devised a scheme to implement general stabilizer quan- 
tum error correction [1] using continuous measurement 
and feedback. Unfortunately an exact implementation 
of the ADL scheme is computationally demanding. For 
an n-qubit code, the procedure requires one to time- 
evolve a 2^-dimensional density matrix for the logical 
qubit alongside the quantum computation [7]. This clas- 
sical information-processing overhead must be performed 
to interpret the continuous-time error syndrome measure- 
ment data and determine how recovery operations, in the 
form of a time-dependent feedback Hamiltonian, should 
be applied. While n is a constant for any particular 
choice of code, even modest codes such as the five-qubit 
code [8, 9] and the seven-qubit Steane code [10] push clas- 
sical computers to their limits. Despite state-of-the art 
experimental capabilities, it would be extremely difficult 
to implement the ADL bit-flip code in practice. Conse- 



quently, Ahn and others have devised alternate feedback 
protocols which are less demanding [11, 12], but perform 
worse than the the original ADL scheme. 

Recently, van Handel and Mabuchi addressed the com- 
putational overhead of continuous-time error syndrome 
detection [13] using techniques from quantum filtering 
theory [14, 15, 16]. They developed an exact, low- 
dimensional model for continuous-time error syndrome 
measurements, but did not go on to treat continuous- 
time recovery. The complication is that any feedback 
Hamiltonian suitable for correcting errors during the syn- 
drome measurements violates the dynamical symmetries 
that were exploited to obtain the low-dimensional filter in 
Ref. [13]. While one might address this complication by 
simply postponing error recovery operations until a point 
where the measurements can be stopped, there may be 
scenarios where it would be preferable to perform error 
recovery in real-time. For example, if the recovery oper- 
ation is not instantaneous, responding to errors as they 
occur might outperform protocols where there are peri- 
ods without any error correction. 

In this paper, we extend the quantum filtering ap- 
proach developed by van Handel and Mabuchi to include 
recovery operations. We consider an error-correcting 
feedback Hamiltonian of the form devised by Ahn, Do- 
herty and Landahl, but our approach readily extends 
to other forms for the feedback. While an exact low- 
dimensional model for continuous-time stabilizer gener- 
ator measurements in the presence of feedback does not 
appear to exist, we devise an approximate filter that is 
still low-dimensional, yet sufficiently accurate such that 
high-quality error correction is possible. 



II. CONTINUOUS-TIME QUANTUM ERROR 
CORRECTION 
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For our purposes, a quantum error correcting code is a 
triple {E, G, R). The quantum operation E : C'^^ ^ C^^ 
encodes k logical qubits in n physical qubits. 5 is a set 



of I = n — k stabilizer generator observables with out- 
comes ±1 that define the error syndrome. R : {±1}^^ \-^ 
£2nx2n ^g ^^^ recovery operation^ which specifies what 
correction should be applied to the physical qubits in re- 
sponse to the syndrome measurement outcomes. 

The particular choice of code {E^Q^R) is usually made 
with consideration for the nature of the decoherence af- 
fecting the physical qubits [17]. For example, the bit- 
flip code (considered by both ADL and van Handel and 
Mabuchi) improves protection against an error channel 
that applies the Pauli a^ operator to single qubits at a 
rate 7. Here, we adopt the notation that Xn represents 
the Pauli ax operator on qubit n, and similarly for Y^ and 
Zn. In the bit-flip code, E encodes k = 1 qubits in n = 3 
qubits by the map a|0) + p\l) ^ a|000) + /3|111). The 
/ = 2 stabilizer generators are gi = ZZI := Gz ® cf^ ^ I 
and g2 = IZZ := I az cr^; each extracts the parity 
of different qubit pairs. The recovery i?, given the out- 
comes of measuring (^1,^2), is defined by (+1, +1) \-^ /, 
(+1, -1) ^ X3, (-1, +1) ^ Xi and (-1, -1) ^ X2. 

In this paper, we focus primarily on the five- qubit- code 
(n = 5, /c = 1) that increases protection against general 
separable channels, and in particular the continuous-time 
symmetric depolarizing channel that applies all three 
Pauli operators to each of the physical qubits at the same 
rate 7. The five-qubit code has / = 4 stabilizer genera- 
tors {XZZXI, IXZZX, XIXZZ, ZXIXZ}. It is also a 
perfect code in that all 16 distinct syndrome outcomes in- 
dicate distinct errors: one corresponding to the no-error 
condition, and one syndrome for each of the three Pauli 
errors on each of the five qubits. We defer to [1, 2] for 
the encoding and recovery procedures for this code. 



A. Continuous-time Stabilizer Generator 
Measurements 

Quantum error correction can be extended to continu- 
ous time by replacing discrete measurements of the sta- 
bilizer generators ^1 , . . . , ^^ with a set of / continuous ob- 
servation processes [7]. This creates i = 1, . . . , / measure- 
ment records 



dQf^ =2^TT[gipt]dt- 



dW, 



(i) 



(1) 



obtained from the encoded qubit with state pt at time 
t. Here, a^: is a constant called the measurement strength 
that depends upon the physical implementation of the 
continuous measurement and the dW^ are independent 
Wiener processes, each with ]E[(iiy^] = and dW^ = dt 
[15, 18]. We do not consider here how one might imple- 
ment the set of / simultaneous stabilizer generator obser- 
vations other than to comment that doing so in an AMO 
technology would likely involve coupling the n physical 
qubits to a set of electromagnetic field modes and then 
performing continuous photodetection on the scattered 
fields. While this model is rather general, we take the 
same hc for each qubit, implying symmetric coupling of 
the qubits. 



By itself, the continuous measurement record is too 
noisy to permit quantum error correction — one must first 
process the measurement data to deal with the presence 
of the noises dW^ . The most straightforward approach 
toward filtering the noise is via an estimate of the full 
logical qubit density operator pt. Generating a full state 
estimation based on the evolving syndrome measurement 
data is accomplished by the quantum filtering equation 
(in its adjoint form with h = 1) [16] 

n I 

dpt = 1Y.T. ^[^^r^]Ptdt + K ^ V[gi]ptdt 

m=l j i=l 

I 

+V^^H[^,]pt [dQ['^ - 2^TT[g,pt]dt 



-i[Ht,pt]dt, 



(2) 



where j G {x, y, z} and the superoperators are defined 
as: V[a]p = apa - p and H[gi]p = gip ^ pgi-2 Tr [gip] p. 
The first term in the filtering equation accounts for the 
action of the continuous-time symmetric depolarizing 
channel while the second accounts for the effect of cou- 
pling the logical qubit to the field used to implement 
the continuous measurements. Conditioning the quan- 
tum state pt on the continuous measurement occurs via 
the third term, which is driven by the innovations pro- 
cesses dQt — 2^/nTr[gipt]dt. The time evolution pt gen- 
erated by a particular noise realization is generally called 
a trajectory. 

The final term in Eq. (2) describes the action of the 
time-dependent feedback Hamiltonian used to implement 
error recovery. Following Ahn, Doherty and Landahl, we 
choose the feedback Hamiltonian to be of the form 






(3) 



Sm) 



which corresponds to applying Pauli operators cr^- to 

each qubit with a controllable strength A^-^ . The policy 

for determining the feedback strengths A^-^^ at each point 
in time should be chosen optimally. Ahn, Doherty, and 
Landahl obtained their feedback policy by defining the 
codespace projector Hq onto the no error states (states 
which are +1 eigenvectors of all stabilizers) and then 
maximizing the codespace fidelity Tr[HoPt]- Assuming 
a maximum feedback strength Amax, the resulting feed- 
back policy is given by setting 



A 



(m) 



AmaxSgn(Tr[-z[Ho,cr]'^V^]) • 



1. Computational Expense 



(4) 



Because this is a closed-loop strategy, the feedback con- 
troller must determine each A^-^ from the evolving mea- 
surement in real time. The utility of feedback in any real 



setting then relies greatly upon the controller's ability to 
integrate the filtering equation rapidly enough to main- 
tain pace with the quantum dynamics of the qubits. For 
the five-qubit code, 1024 — 1 real parameters are needed 
to represent the density matrix. We found that stable 
numerical integration [19] of even a single trajectory re- 
quired approximately 36 seconds on a 2.1 GHz desktop 
computer {jdt ^ 10~^ over a timespan [0, O.257]). This is 
far from adequate for use in an actual feedback controller 
even in state-of-the-art experiments. 

Moreover, Eq. (2) is a nonlinear filter, and for such 
filters it is rarely possible to evaluate even qualitative 
properties analytically. One must then average over an 
appreciable number of trajectories to find the expected 
behavior of quantities such as the codespace fidelity as a 
function of time. For the five-qubit code, our integrator 
requires approximately 10 hours to simulate 1000 trajec- 
tories. 



B. Reduced-Dimensional Filters 

Considering that the syndrome measurements yield in- 
formation about correlations between qubits and not in- 
formation about the individual states of the qubits, one 
can imagine that propagating the full density matrix 
is excessive. Indeed, the ADL scheme only makes use 
of the projection of pt onto the codespace, generating 
the same feedback policy regardless of which state po in 
the codespace is initially chosen. It is reasonable to ex- 
pect that a lower dimensional model could track solely 
the information extracted from the syndrome measure- 
ments. This is exactly the premise used by van Han- 
del and Mabuchi to obtain a low-dimensional model of 
continuous-time stabilizer generator measurements (in 
the absence of feedback) [13]. They formulate the prob- 
lem as a graph whose vertices correspond to syndromes 
and whose edges reflect the action of the error model. 
The filtering problem is then reduced to tracking the 
node probabilities, i.e., the likelihoods for the qubit to 
be described by each of the various syndrome conditions. 
Dynamical transitions occur between the syndromes due 
to the error channel, and the filter works to discern these 
transitions from the stabilizer measurement data. 

For an (E^Q^R) code, van Handel and Mabuchi define 
a set of projectors onto the distinct syndrome spaces. For 
the five-qubit code, there are 16 such projectors; Hq is the 



(m) 



af^Hoa^^^) 



codespace projector as before and H 

are projectors onto states with a syndrome consistent 

with a aj error on qubit m. Forming the probabilities 



(m) 



^Tr[n^-) 



Pt 



S^) 



(5) 



into a vector p^ and computing dp- ^ from the full dy- 
namics leads to the reduced filter 

I 
dY>t = Apt dt + 2^/^Y^{Hi - hi^p, /)p, dWt (6) 



with Krs = 7(1 — lQ5rs)^ hj'^ the outcome of measuring 
gi on H^."^^ and Hi = diagh^ (Eq. (4) in Ref. [13]). The 

equations for pj^^ are closed and encapsulate all the in- 
formation that is gathered from measuring the stabilizer 
generators. Equation (6) is an example of a Wonham fil- 
ter^ which is the classical optimal filter for a continuous- 
time finite-state Markov chain with an observation pro- 
cess driven by white noise [20] . Further discussion of the 
Wonham filter and its use in conjunction with discrete- 
time error correction can be found e.g., in Ref. [13]. 



III. CONTINUOUS-TIME QUANTUM 
FILTERING WITH FEEDBACK 



We now extend Eq. (6) to include a feedback Hamil- 
tonian suitable for error recovery. Following van Handel 
and Mabuchi's lead, we see that Eq. (6) was derived by 

taking dp^ = Tr[H^-^^(ipt] for a basis which closed un- 
der the dynamics of the continuous syndrome measure- 
ment. One hope is that simply adding the feedback term 
in by calculating Tr[— iA^^^H^-^^[(j^^ , p^]] also results in 
a set of closed equations. However, that is not the case 
when using the basis of the sixteen syndrome space pro- 
jectors H^^K Specifically, [H^"^\cr^^^] cannot be written 
as a linear combination of syndrome space projectors. 
This is not surprising as the feedback Hamiltonian term 
under consideration is the only term which generates uni- 
tary dynamics. 

Inspired by the form of the commutator between the 
feedback and the syndrome space projectors, we define 
feedback coefficient operators 



H 



(m) 



(+z or + l)cr^ 



5n(-)^05^ 



(7) 



where c is an arbitrarily chosen index used to distinguish 
the z or 1 prefactor and combination of Pauli matrices 
which sandwich the syndrome space projector H^-^K For 
the five-qubit code, the syndrome projectors are simply 
those operators which have the 1 prefactor and 10 iden- 
tity matrices. The corresponding feedback coefficient is 

(m) • ■ 



P 



Tr [H^-^^pt] • If we then iterate the dynamics of the 



Sm) 



filter (2) by calculating pj^J starting from the syndrome 
space projectors, we find that each feedback Hamilto- 
nian term generates pairs of feedback coefficient terms. 
For example, calculating the dynamics due to feedback 
Xi on Hq generates two feedback coefficient operators: 
Ho,o = ^HqXi and Ho,i = iXiUo. We must then deter- 
mine the dynamics for these first level feedback coeffi- 
cients. This will include calculating the Y^ feedback on 
Ho,i, which generates second level feedback coefficients 
Ho,2 = -^i^sHo and Ho,3 = XiUqY^. Continuing to it- 
erate feedback coefficient terms, we find that an addi- 
tional 1008 distinct pj^^ terms are needed to close the 
dynamics and form a complete basis. Adding in the ini- 
tial 16 syndrome space projectors gives a 1024 dimen- 
sional basis — clearly no better than propagating the full 
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FIG. 1: Non-zero matrix elements of (a) untruncated and (b) truncated filter. Blue squares correspond to decoherence terms, 
red crosses correspond to measurement terms and green dots correspond to feedback terms. Note the difference in dimension 
of the matrices. (Color online.) 



density matrix. However, it is now relatively easy to cal- 
culate the feedback strengths, which depend only on pairs 
of first-level feedback coefficients. For example, from Eq. 



.(1) 



(4) we find that A^^ = A^ax sgn (-po,o +Po,i), where 
Po,o = Tr[no,oP^] and po,i = Tr[no,iPt] are first-level 
coefficients developed earlier in the paragraph. 



A. Approximate Filter for the Five-Qubit Code 

Although the dimension of the alternate basis is no 
smaller than the dimension of the full density matrix, the 
structure of the filter represented in the alternate basis 
provides a manner for interpreting the relative impor- 
tance of the p^ feedback coefficients. This is best seen 
graphically in Fig. 1(a), which superimposes the non- 
zero matrix elements coming from the noise, measure- 
ment and feedback terms. Both measurement and noise 
are block diagonal as expected; it is the feedback that 
couples blocks together in a hierarchical fashion. This hi- 
erarchy can be parameterized by the number of "feedback 
transitions" which connect a given feedback coefficient to 
the syndrome space block. For example, the upper left 
block, which corresponds to the syndrome space projec- 
tors, is connected via feedback terms to the first level 
feedback block, whose feedback coefficients are each one 
feedback transition away from the syndrome space block. 
In turn, the first level block is then connected to a second 
level feedback block, whose feedback coefficients are two 
feedback transitions away from the syndrome block. 



Given that the initial state starts within the codespace 
and given that feedback is always on, the feedback coef- 
ficients that are more than one feedback transition from 
the syndrome space block should be vanishingly small. 
Limiting consideration to these first two blocks, we also 
find that pairs of feedback coefficients couple identically 
to the syndrome space block. For example, we find that 
—iXiIlo and HIqXi couple to syndrome space projec- 
tors identically. This is not surprising, as these two 
terms comprise the commutator that results from the 
Xi feedback Hamiltonian. However, outside the first 
level of feedback transitions, the matrix elements of these 
feedback coefficients differ. Additionally, feedback coeffi- 
cients involving feedback Hamiltonians which correspond 
to a syndrome error on the codespace projector are re- 
lated as 






iUna["''> = -m 



(to) (to) 



■"'j 



{m)jj{m) 



n 



• (8) 



For the feedback coefficient examples just mentioned, this 
relation is -zXiH^^^+zH^^^Xi = -zHo^i+^^ino- Trun- 
cating the dynamics to include only the first level of feed- 
back and combining distinct feedback coefficients which 
act identically within this block results in the matrix of 
Fig. 1(b) over only 136 basis elements. Note that the con- 
troller now only needs to reference a single basis element 
for calculating a given feedback strength A^-^ . 



B. Approximate Filter for General Codes 

Our truncation scheme generalizes for reducing the 
dimensionality of the quantum filter for an arbitrary 
(E^Q^R) code. Such a filter for an |n, kj quantum error- 
correcting code [2] has the same form as Eq. (2), but 
involves n physical qubits and I = n — k continuous- 
time stabilizer generator measurements. In the following, 
we assume the continuous-time symmetric depolarizing 
channel, though it should be straightforward to extend 
to other noise models. For a non-perfect, non-degenerate 
code, there are a total of 2"^"^ stabilizer generator mea- 
surement outcomes, but only 3n + 1 will be observed for 
the given noise channel. For a perfect, non-degenerate 
code (2'^~^ = 3n + l), all possible syndrome outcomes are 
observed. In either case, given the observable syndrome 
outcomes, we can define 3n-\-l syndrome space projectors 
and 3n feedback parameters needed for recovery. Degen- 
erate codes require fewer than 3n recovery operations, as 
distinct actions of the noise channel give rise to identical 
errors and recovery operations. The degeneracy depends 
greatly on the particular code, so we merely note that de- 
generate codes will require fewer syndrome space projec- 
tors and feedback parameters than their non-degenerate 
relatives. 

Once we determine the syndrome space projectors and 
feedback parameters for the code, we can introduce feed- 
back coefficient operators of the form of (7) but over n 
qubits. A truncated filter is constructed as follows. 

1. Close the dynamics of the 3n + 1 syndrome space 
projectors by introducing 6n(3n+l) first-level feed- 
back terms (2 feedback coefficients per commutator 
in each of the 3n feedback Hamiltonians). 

2. Close the dynamics of the first-level feedback terms, 
truncated to a basis of syndrome space and first- 
level feedback terms, i.e. ignore any second-level 
terms which were not defined in step 1. The proce- 
dure involves 0{n^) steps, since one must at least 
examine each of the 6n feedback commutator terms 
for all first-level feedback coefficients. 

3. Each of the 3n + 1 syndrome space projectors is 
related to pairs of first-level feedback coefficients 
that comprise a commutator. There is an addi- 
tional factor of degeneracy between syndrome space 
projectors and feedback coefficients which involve 
the same Pauli matrix [c.f., Eq. (8)]. A similarity 
transform is used to eliminate these redundancies, 
leaving (3n + l) + (3n + l)6n/4 = ^ (2 + 9n(n + 1)) 
basis elements in the fully truncated filter. 

The truncated filter requires only 0{ii?) basis elements, 
as compared to the 4^ parameters for the full density ma- 
trix. Additionally, the feedback strengths in Eq. (4) are 
readily calculated from the combined first-level feedback 
coefficients. The truncation process is depicted schemat- 
ically in the left half of Fig. 2. The right half of the figure 



Retained Elements 



Syndrome Space 
3n Feedback Hamiltonians 




3n Feedback [Hamiltonians 



^^2^ 



Level Feedback Block 



lock^ 



Schennatic Overview 



Retained Terms 



^^<^^>;;d&^ 



1st Level Feedback 
(Xi Terms) 



2"d Level Feedback 
(75 Terms) 

Exannple Terms 



FIG. 2: On the left, a schematic diagram of truncating the 
filter to only syndrome space and first level feedback blocks. 
On the right, just a few of the 1024 feedback coefficients of 
the five-qubit code representing the different feedback block 
levels. 



gives exaniples of a few of the 1024 terras involved in the 
truncation procedure for the five-qubit code. 



C. Numerical Simulation 

Since the truncated filter is also nonlinear, it is difficult 
to provide analytic bounds on possible degradation in 
performance. However, we can easily compare numerical 
simulation between feedback controllers which use the 
full or truncated filter. In fact, the dynamics should be 
close for the same noise realizations, indicating that they 
should be close per trajectory. 

In order to analyze the feedback controller's perfor- 
mance, the full filter Eq. (2) is used to represent the 
underlying physical system. The feedback controller was 
modeled by simultaneously integrating the truncated fil- 
ter, driven by the measurement current from the full fil- 
ter. The feedback controller then calculated the feedback 
strengths which were fed back into the full filter. The dy- 
namics described by the full filter were then used to com- 
pute the codespace fidelity. Using a predictor-corrector 
stochastic differential equation (SDE) integrator [19] and 
varying n and Amax over a wide range, we found essen- 
tially indistinguishable performance between the full and 
truncated filters. Using n = IOO7 and Amax = 2OO7 as 
representative parameters. Figure 3(a) demonstrates this 
general behavior by comparing the average codespace fi- 
delity of a handful of trajectories using the different fil- 
ters. Integrating an individual trajectory takes approx- 
imately 39.5 seconds using a 2.1 GHz PowerPC proces- 
sor. Integrating the full filter alone takes approximately 
36 seconds, while integrating the truncated filter alone 
takes approximately 3.5 seconds. 

In addition to showing the identical performance of the 
full and truncated filters. Fig. 3(a) also shows the loss in 
performance if one were to truncate further. The 31 di- 
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FIG. 3: Numerical simulations of the five qubit code to assess the average code space fidelity. Plot (a) compares the codespace 
fidelity (averaged over 10 trajectories) for filters with different levels of truncation: the full (1024-dimensional) and first-level 
truncated (136-dimensional) filters are essentially identical. Plot (b) shows the codespace fidelity averaged over 2,000 trajectories 
using the truncated 136-dimensional filter for error correction. (Simulation parameters: Amax = 2OO7 and k = IOO7.) (Color 
online.) 



mensional filter is comprised of the 16 syndrome projec- 
tors and the 15 feedback coefficients which have non-zero 
feedback matrix elements with the codespace IIo. These 
are the only elements explicitly needed to calculate the 
feedback strengths in Eq. 4. This filter fails because it 
tacitly assumes the action of feedback on the codespace is 
more "important" than on the other 15 syndrome spaces. 
Since feedback impacts all syndrome spaces equally, we 
need to retain those terms in order to properly maintain 
syndrome space probabilities. Intuitively, this suggests 
that the 136 dimensional filter is the best we can do us- 
ing this heuristic truncation strategy. For reference. Fig. 
3(b) shows the average codespace fidelity of 2000 trajec- 
tories when using the truncated filter. 



1. Comparison with Discrete Error Correction 

Given the success of the truncation scheme, we now 
compare the performance of feedback- assisted error cor- 
rection to that of discrete-time error correction for the 
five-qubit code. The discrete model considers qubits ex- 
posed to the depolarizing channel 



dp, 



'discrete 



n=5 

E 

j=x,y,z m=l 



E 



^[^.-"^VdiscreteC^^ (9) 



up to a time t, after which discrete- time error correction 
is performed. The solution of this master equation can 
be explicitly calculated using the ansatz 

5 

Pdiscrete{t) = Yl ^ ae{t)PpoP, (10) 

e=0 P;pw{P)=e 

where P is a tensor product of Pauli matrices and the 
identity. The function pw{P) gives the Pauli weight 



of a matrix, defined as the number of ax^Cy^ and a^ 
terms in the tensor representation. Thus, ao{t) is the 
coefficient of po and similarly ai{t) is the coefficient 
of all single qubit errors from the initial state, e.g., 
XIIII{po)XIIII, IIZII{po)IIZIL 

The codespace fidelity considered earlier is not a useful 
metric for comparison, as discrete-time error correction is 
guaranteed to restore the state to the codespace. Follow- 
ing Ahn, Doherty and Landahl, we instead use the code- 
word fidelity Fcw{t) := Tr[pop(^)], which is a measure 
relevant for a quantum memory. Since error correction is 
independent of the encoded state, we choose the encoded 
|0) state as a fiducial initial state. Given that the five- 
qubit code protects against only single qubit errors, we 
find that after error correction at time t, the codeword 
fidelity for discrete-time error correction is 



TT'discrete ^ /j.\ 



256 



-ai{t) 



which asymptotes to 1/64. This limit arises because 
prior to the stabilizer generator measurements, the noise 
pushes the state to the maximally mixed state, which 
is predominately composed of the a2{t) through a^{t) 
terms. 

The feedback codeword fidelity i^^Jeedback ^^^ calcu- 
lated by integrating both the full quantum filter (2), rep- 
resenting the underlying system of qubits, and the trun- 
cated filter, representing the feedback controller. Again, 
we chose hz = IOO7, Amax = 2OO7 and dt = 10~^7 and 
used the same SDE integrator described above. Fig- 
ure 4 shows the average of i^-f^edback ^^^^ 2000 tra- 
jectories, demonstrating that there are regimes where 
feedback-assisted error correction can significantly out- 
perform discrete-time error correction. Feedback- assisted 
error correction appears to approach an asymptotic code- 
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FIG. 4: Comparison between continuous-time and discrete- 
time error correction for the five-qubit code. For the 
continuous-time error correction simulations, the codeword 
fidehty was averaged over 2,000 trajectories with n = IOO7 
and Amax = 2OO7. (Color online.) 



word fidelity greater than what would be obtained by de- 
coherence followed by discrete-time error correction. Due 
to the nonlinear feedback, we were unable to calculate an 
analytic asymptotic expression for the continuous-time 
strategy. Nonetheless, the improved performance for the 
timespan considered suggests that better quantum mem- 
ory is possible using the feedback scheme. 



for the original Ahn, Doherty and Landahl procedure [7] . 
By numerical simulation of the five-qubit code, we have 
demonstrated the viability of such a filter for a quantum 
memory protecting against a depolarizing noise channel. 
Moreover, in all our simulations, this performance is in- 
distinguishable from that of the computationally more 
demanding filter of the ADL style. 

In systems where recovery operations are not instan- 
taneous relative to decoherence, consideration suggests 
that it is desirable to perform syndrome measurement, 
recovery, and logic gates simultaneously. However, it 
is not immediately clear how gates impact the feedback 
controller. Indeed, if a Hamiltonian is in the code's nor- 
malizer, the continuous-time feedback protocol and its 
performance are unchanged. Though a universal set of 
such Hamiltonians can be found, it might be desirable to 
find universal gates which have physically simple inter- 
actions. Future work involves finding such gate sets and 
developing a framework for universal quantum comput- 
ing. Additional issues of fault-tolerance and robustness 
could then be explored within such a universal setup. Ex- 
ploring feedback error correction in the context of spe- 
cific physical models will provide opportunities to tailor 
feedback strategies to available control parameters and 
salient noise channels. Such systems might allow the cal- 
culation of globally optimal feedback control strategies. 
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